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A review of the development of the concept of dark matter is 
given. The dark matter story passed through several stages on 
its way from a minor observational puzzle to a major challenge 
for theory of elementary particles. 



I begin the review with the description of the discovery of 
the mass paradox in our Galaxy and in clusters of galaxies. First 
hints of the problem appeared already in 1930s and later more 
observational arguments were brought up, but the issue of the 
mass paradox was mostly ignored by the astronomical commu- 
nity as a whole. In mid 1970s the amount of observational data 
was sufficient to suggest the presence of a massive and invisi- 
ble population around galaxies and in clusters of galaxies. The 
nature of the dark population was not clear at that time, but the 
hypotheses of stellar as well as of gaseous nature of the new pop- 
ulation had serious difficulties. These difficulties disappeared 
when non-baryonic nature of dark matter was suggested in early 
1980s. 

The final break through came in recent years. The systematic 
progress in the studies of the structure of the galaxies, the studies 
of the large scale structure based on galaxy surveys, the analysis 
of the structure formation after Big Bang, the chemical evolu- 
tion of the Universe including the primordial nucleosynthesis, 
as well as observations of the microwave background showed 
practically beyond any doubt that the Universe actually contains 
more dark matter than baryonic matter! In addition to the pres- 
ence of Dark Matter, recent observations suggest the presence 
of Dark Energy, which together with Dark Matter and ordinary 
baryonic matter makes the total matter/energy density of the 
Universe equal to the critical cosmological density. Both Dark 
Matter and Dark Energy are the greatest challenges for modern 
physics since their nature is unknown. 

There are various hypothesis as for the nature of the dark 
matter particles, and generally some form of weakly interactive 
massive particles (WIMPs) are strongly favored. These parti- 
cles would form a relatively cold medium thus named Cold Dark 
Matter (CDM). The realization that we do not know the nature 
of basic constituents of the Universe is a scientific revolution 
difficult to comprehend, and the plan to hunt for the dark matter 
particles is one of the most fascinating challenges for the future. 

1 Dark Matter problem as a scientific 
revolution 

Almost all information on celestial bodies comes to us via pho- 
tons. Most objects are observed because they emit light. In other 
cases, like for example in some nebulae, we notice dark regions 
against otherwise luminous background which are due to ab- 
sorption of light. Thus both light absorption and light emission 
allow to trace the matter in the Universe, and the study goes 
nowadays well beyond the optical light. Modern instruments 
have first detected photon emission from astronomical bodies in 
the radio and infrared regions of the spectrum, and later also in 
the X-ray and gamma-ray band, with the use of detectors in- 
stalled in space. 



Presently available data indicate that astronomical bodies of 
different nature emit (or absorb) photons in very different ways, 
and with very different efficiency. At the one end there are ex- 
tremely luminous supernovae, when a single star emits more en- 
ergy than all other stars of the galaxy it belongs to, taken to- 
gether. At the other extreme there are planetary bodies with a 
very low light emission per mass unit. The effectiveness of the 
emissivity can be conveniently described by the mass-to-light 
ratio of the object, usually expressed in Solar units in a fixed 
photometric system, say in blue (B) light. The examples above 
show that the mass-to-light ratio Mj L varies in very broad range. 
Thus a natural question arises: Do all astronomical bodies emit 
or absorb light? Observations carried out in the past century 
have led us to the conclusion that the answer is probably NO. 

Astronomers frequently determine the mass by studying the 
object emission. However, the masses of astronomical bodies 
can be also determined directly, using motions of other bodies 
(considered as test particles) around or within the body under 
study. In many cases such direct total mass estimates exceed the 
estimated luminous masses of known astronomical bodies by a 
large fraction. It is customary to call the hypothetical matter, 
responsible for such mass discrepancy, Dark Matter. 

The realization that the presence of dark matter is a serious 
problem which faces both modern astronomy and physics grew 
slowly but steadily. Early hints did not call much attention. 

The first indication for the possible presence of dark mat- 
ter came from the d ynamical stu dy of our Galaxy. Dutch as- 
tronomer Jan Henrik lOorJ lll932h analysed vertical motions of 
stars near the plane of the Galaxy and calculated from these data 
the vertical acceleration of matter. He also calculated the vertical 
acceleration due to all known stars near the Galactic plane. His 
result was alarming: the density due to known stars is not suf- 
ficient to explain vertical motions of stars - there must be some 
unseen matter near the Galactic plane. 

The second observation was made by Fritz IZwicky (1933). 
He measured radial velocities of galaxies in the Coma cluster of 
galaxies, and calculated the mean random velocities in respect 
to the mean velocity of the cluster. Galaxies move in clusters 
along their orbits; the orbital velocities are balanced by the total 
gravity of the cluster, similar to the orbital velocities of planets 
moving around the Sun in its gravitation field. To his surprise 
Zwicky found that orbital velocities are almost a factor of ten 
larger than expected from the summed mass of all galaxies be- 
longing to the cluster. Zwicky concluded that, in order to hold 
galaxies together in the cluster, the cluster must contain huge 
amounts of some Dark (invisible) matter. 

The next hint of the dark matter existence came from cos- 
mology. 

One of the cornerstones of the modern cosmology is the con- 
cept of an expanding Universe. From the expansion speed it is 
possible to calculate the critical density of the Universe. If the 
mean density is less than the critical one, then the expansion con- 
tinues forever; if the mean density is larger than the critical, then 
after some time the expansion stops and thereafter the Universe 
starts to collapse. The mean density of the Universe can be esti- 
mated using masses of galaxies and of the gas between galaxies. 
These estimates show that the mean density of luminous matter 
(mostly stars in galaxies and interstellar or intergalactic gas) is 
a few per cent of the critical density. This estimate is consistent 
with the constraints from the primordial nucleosynthesis of the 
light elements. 

Another cornerstone of the classical cosmological model is 



the smooth distribution of galaxies in space. There exist clusters 
of galaxies, but they contain only about one tenth of all galax- 
ies. The majority of galaxies are more or less randomly dis- 
tributed and are called field galaxies. This conclusion is based 
on counts of galaxies at various magnitudes and on the distribu- 
tion of galaxies in the sky. 

Almost all astronomical data fitted well to these classical 
cosmological paradigms until 1970s. Then two important anal- 
yses were made which did not match the classical picture. In 
mid 1970s first redshift data covering all bright galaxies were 
available. These data demonstrated that galaxies are not dis- 
tributed randomly as suggested by earlier data, but form chains 
or filaments, and that the space between filaments is practically 
devoid of galaxies. Voids have diameters up to several tens of 
megaparsecs. 

At this time it was already clear that structures in the Uni- 
verse form by gravitational clustering, started from initially 
small fluctuations of the density of matter. Matter "falls" to 
places where the density is above the average, and "flows away" 
from regions where the density is below the average. This grav- 
itational clustering is a very slow process. In order to form 
presently observed structures, the amplitude of density fluctu- 
ations must be at least one thousandth of the density itself at the 
time of recombination, when the Universe started to be transpar- 
ent. The emission coming from this epoch was first detected in 
1965 as a uniform cosmic microwave background. When finally 
the fluctuations of this background were measured by COBE 
satellite they appeared to be two orders of magnitude lower than 
expected from the density evolution of the luminous mass. 

The solution of the problem was suggested independently 
by several theorists. In early 1980s the presence of dark matter 
was confirmed by many independent sources: the dynamics of 
the galaxies and stars in the galaxies, the mass determinations 
based on gravitational lensing, and X-ray studies of clusters of 
galaxies. If we suppose that the dominating population of the 
Universe - Dark Matter - is not made of ordinary matter but of 
some sort of non-baryonic matter, then density fluctuations can 
start to grow much earlier, and have at the time of recombination 
the amplitudes needed to form structures. The interaction of 
non-baryonic matter with radiation is much weaker than that of 
ordinary matter, and radiation pressure does not slow the early 
growth of fluctuations. 

The first suggestion for the non-baryonic matter were parti- 
cles well known at that time to physicists - neutrinos. However, 
this scenario soon lead to major problems. Neutrinos move with 
very high velocities which prevents the formation of small struc- 
tures as galaxies. Thus some other hypothetical non-baryonic 
particles were suggested, such as axions. The essential property 
of these particles is that they have much lower velocities. Be- 
cause of this the new version of Dark Matter was called Cold, 
in contrast to neutrino-dominated Hot Dark Matter. Numeri- 
cal simulations of the evolution of the structure of the Universe 
confirmed the formation of filamentary superclusters and voids 
in the Cold Dark Matter dominated Universe. 

The suggestion of the Cold Dark Matter has solved most 
problems of the new cosmological paradigm. The actual nature 
of the CDM particles is still unknown. Physicists have attempted 
to discover particles which have properties needed to explain the 
structure of the Universe, but so far without success. 

One unsolved problem remained. Estimates of the matter 
density (ordinary + dark matter) yield values of about 0.3 of the 
critical density. This value - not far from unity but definitely 



smaller than unity - is neither favoured by theorists nor by the 
data, including the measurements of the microwave background, 
the galaxy dynamics and the expansion rate of the Universe ob- 
tained from the study of supernovae. To fill the matter/energy 
density gap between unity and the observed matter density it 
was assumed that some sort of vacuum energy exists. This as- 
sumption is not new: already Einstein added to his cosmological 
equations a term called the Lambda-term. About ten years ago 
first direct evidence was found for the existence of the vacuum 
energy, presently called Dark Energy. This discovery has filled 
the last gap in the modern cosmological paradigm. 

In the International Astronomical Union (IAU) s ymposium 
on Dark Matter in 1985 in Princeton. iTremaine ( 1987 ) charac- 
terised the discovery of the dark matter as a typic al scientific r ev- 
olution, connected with changes of paradigms. iKuhnl dl970h in 
his book The Structure of Scientific Revolutions discussed in de- 
tail the character of scientific revolutions and paradigm changes. 
There are not so many areas in modern astronomy where the de- 
velopment of ideas can be described in these terms, thus we shall 
discuss the Dark Matter problem also from this point of view. 
Excellen t reviews on the dark matt e r and related problems are 
iven b v lFaber & Gallagher! d 1979b iTrimbld d 1987b. ISrednickH 
1990b. [Turnerl dl99lb ISilld dl992b Ivan den Berghl d200l | ) . 
Ostriker & Steinhardtl 12003b. Reesl d2003h . iTurnerj (12003b . 

d2008b. see also proceed- 



(1987) 



Tegmark et al.l d2006b. iFrieman et alj _ 

ings b y Longair & Einastol ~( 1978 ). and Kormendv & Knappl 



2 Early evidence of the existence of dark 
matter 

2.1 Local Dark Matter 

The dynamical density of matter in the Solar vicinity can be es- 
timated using vertical oscillations of stars around the galactic 
plane. The orbital motions of stars around the galactic center 
play a muc h smaller role in determining the local density. Ernst 
Opik d 1 9 1 5D found that the summed contribution of all known 
stellar populations (and interstellar gas) is sufficient to explain 
the vertical oscillations of stars - in other words, there is no 
need to assume the e xistence of a dar k population. A similar 
analysis was made by iKaptevnl d 1922b . who first used the term 
"Dark Matter" to denote invisible matter which existence is sug- 
gested by its gravity only. Both Opik and Kapteyn found that the 
amount of invisible matter in the Solar neighbourho od is s mall. 

Another conclusion was obtained by Jan lOortl dl932l) . His 
analysis indicated that the total density, found from dynamical 
data, exceeds the density of visible stellar populations by a factor 
of up to 2. This limit is often called the Oort limit. This result 
means that the amount of invisible matter in the Solar vicinity 
should be approximately equal to the amount of visible matter. 

The local density of matt er has been redetermin ed by var- 
ious authors man y time s. Kuzmin dl952bl 1 1955b and his 
stude nts Heino lEelsalul d 1959b and Mihkel IJoeveerl dl972 , 
1974) confirmed the earlier result by Opik. A number 



of other astronom e rs, in cl uding m o re rec ently lOoril (I 1960b 



Bahcall & Soneiral dl980b ; iBahcalll {[984), found results 



in 



agreement with Oorts original result. This discussion was open 
until recently; we will describe the present conclusions below. 

For long time no distinction between local and global dark 
matter was made. The realisation, that these two types of dark 



matter have very different properties and nature came from the 
detailed study of galactic models, as we shall discuss below. 

2.2 Global Dark Matter - clusters, groups and 
galaxies 



A different mass discrepancy was found by Fritz IZwick v ( 1933). 
He measured redshifts of galaxies in the Coma cluster and found 
that the velocities of individual galaxies with respect to the clus- 
ter mean velocity are much larger than those expected from the 
estimated total mass of the cluster, calculated from masses of 
individual galaxies. The only way to hold the cluster from rapid 
expansion is to assume that the cluster contains huge quantities 
of some invisible dark matter. According to his estimate the 
amount of dark matter in this cluster exceeds the total mass of 
cluster galaxies at least tenfolds, probably even more. As char- 
acteristic in scientific revolutions, early indications of problems 
in current paradigms are ignored by the community, this hap- 
pened also with the Zwicky's discovery. 

The next step in the study of masses of systems of galax- 
ies was made by iKahn & Woltierl dl959l) . They paid attention 
to the fact that most galaxies have positive redshifts as a result 
of the expansion of the Universe; only the Andromeda galaxy 
(M31) has a negative redshift of about 120 km/s, directed to- 
ward our Galaxy. This fact can be explained, if both galaxies, 
M31 and our Galaxy, form a physical system. A negative ra- 
dial velocity indicates that these galaxies have already passed 
the apogalacticon of their relative orbit and are presently ap- 
proaching each other. From the approaching velocity, the mutual 
distance, and the time since passing the perigalacticon (taken 
equal to the present age of the Universe), the authors calcu- 
lated the total mass of the double system. They found that 
M tot > 1.8 x 10 12 M Q . The conventional masses of the Galaxy 
and M31 are of the order of 2 x 10 11 M . In other words, the 
authors found evidence for the presence of additional mass in 
the Local Group of galaxies. The authors suggested that the 
extra mass is probably in the form of hot gas of temperature 
about 5x1 5 K. Using more modern data lEinasto & Lvnden-Bell 
dl982b made a new estimate of the total mass of the Local 
Group, using the same approach, and found the total mass of 
4.5 + 0.5 x 1O 12 M . This estimate is in good agreement with 
new determinations of the sum of masses of M3 1 and the Galaxy 
including their dark halos (see below). 

A certain discrepancy was also detected between masses of 
individual galaxies and masses of pairs and groups of galaxies. 
The conventional approach for the mass determination of pairs 
and groups of galaxies is statistical. The method is based on the 
virial theorem and is almost identical to the procedure used to 
calculate masses of clusters of galaxies. Instead of a single pair 
or group often a synthetic group is used consisting of a number 
of individual pairs or groups. These determinations yield for the 
mass-to-light ratio (in blue light) the values M/Lb = 1 ... 20 for 
spiral galaxy dominated pairs, and M/L B = 5 ... 90 for e llipti- 
cal ga laxy dominated pairs (for a review see iFaber & Galla gher 
dl979b ). These ratios are larger than found from local mass indi- 
cators of galaxies (velocity dispersions at the center and rotation 
curves of spiral galaxies). However, it was not clear how seri- 
ous is the discrepancy between the masses found using global or 
local mass indicators. 
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Figure 1 : The rotation curve of M3 1 bv lRoberts & Whitehursti(ll975b . The filled triangles show the optical data from lRubin & Ford 
(1970), the filled circles show the 2 1 -cm measurements made with the 300-ft radio telescope (reproduced by permission of the AAS 
and the author). 



2.3 Rotation curves of galaxies 

Another problem with the distribution o f mass and ma ss-to-light 
ratio was detected in spiral galaxies. iBabcockl ( 1939b obtained 
spectra of the Andromeda galaxy M31, and found that in the 
outer regions the galaxy is rotating with an unexpectedly high 
velocity, far above the expected Keplerian velocity. He inter- 
preted this result either as a high mass -to-light rati o in the pe- 
riphery or as a strong dust absorption. lOortl (119401) studied the 



rotation and surface brightness of the edge-on SO galaxy NGC 
3115, and foun d in the outer region s a mass- to-light ratio ~ 250 . 
Subsequently. iRubin & Fordl d 1 97Clb and lRoberts & Rotsl d 1973b 
extended the rotation curve of M3 1 up to a distance ~ 30 kpc, us- 
ing optical and radio data, respectively. The rotation speed rises 
slowly with increasing distance from the center of the galaxy 
and remains almost constant over radial distances of 16-30 kpc, 
see Fig.Q] 

The rotation data allow to determine the distribution of mass, 
and the photometric data - the distribution of light. Compar- 
ing both distributions one can calculate the local value of the 
mass-to-light ratio. In the periphery of M3 1 and other galaxies 
studied the local value of M/L, calculated from the rotation and 
photometric data, increases very rapidly outwards, if the mass 
distribution is calculated directly from the rotation velocity. In 
the periphery old metal-poor halo-type stellar populations dom- 
inate. These metal -poor populations have a low M/L » 1 (this 
value can be checked directly in globular clusters which con- 
tain similar old metal-poor stars as the halo). In the peripheral 
region the luminosity of a galaxy drops rather rapidly, thus the 
expected circular velocity should decrease according to the Ke- 
plerian law. In contrast, in the periphery the rotation speeds of 
galaxies are almost constant, which leads to very high local val- 
ues of M/L > 200 near the last points with a measured rotational 
velocity. 

Two possibilities were suggested to solve this controversy. 
One possibility is to identify the observed rotation velocity with 
the circular velocity, but this leads to the presence in galaxies 
of an extended population with a very high M/L. The other 
possibility is to assume that in the periphery of galaxies there 



exist non-circular motions which distort the rotation velocity. 

To make a choice between the two possibilities for solving 
the mass discrepancy in galaxies more detailed models of galax- 
ies were needed. In particular, it was necessary to take into ac- 
count the presence in galaxies stellar populations with different 
physical properties (age, metal content, colour, M/L value). 

2.4 Mass paradox in galaxies from Galactic mod- 
els 

Classical models of elliptical galaxies were found from luminos- 
ity profiles and calibrated using either central velocity disper- 
sions, or motions of companion galaxies. The luminosity pro- 
files of disks were often approximated by an e xponential law, 
and b ulge and halo dominated ellipticals by the Ide Vaucouleurs 
dl953bh law. 

Models of spiral galaxies were constructed using rotation ve- 
locities. As a rule, the rotation ve locity was approx imated by 
some simple formula, such as the iBottlingerl (11933b law, or a 
polynomial. The other possibility was to approximate the spa- 
tial density (calculated from t he rotation data ) by a sum of ellip- 
soids of constant density (the Schmidt (1956) model). In the first 
case there exists a danger that, if the velocity law is not chosen 
well, then the density in the periphery of the galaxy may have 
unrealistic values (negative density or too high density, leading 
to an infinite total mass). If the model is built by superposition 
of ellipsoids of constant density, then the density is not a smooth 
function of the d istance from the center of the galaxy. To avoid 
these difficulties iKuzminl i 1 952al 1 1 956) developed models with 
a continuous change of the spatial density, and applied the new 
technique to M3 1 and our Galaxy. His method allows to apply 
this approach also for galaxies consisting of several populations. 

A natural generalisation of classical galactic models is the 
use of all available observational data for spiral and elliptical 
galaxies, both photometric data on the distribution of colour and 
light, and kinematic al data on the rotation and/or velocity disper- 
sion. Further, it is natural to apply identical methods for mod- 
eling of galaxies of different morphological type (including our 



own Galaxy), and to describe explicitly all major stellar popu- 
lations, such as the bulge, the disk, the halo, as well as the flat 
population in spiral galaxies , consisting of young stars and in- 



terstellar gas (lEinastoU l965). 



Multi-component models for spiral an d elliptical gala xies 
using photometric data were constructed by iFreemanl (1 1 970h . To 
combine photometric and kinematic data, mass-to-light ratios of 
galactic populations are needed. Luminosities and colours of 
galaxies in various photometric systems result from the physical 
evolution of stellar populations that can be m odeled. The ch emi- 
cal evolution of galaxies w as investigated by rrTiislev ( 19681) and 
ICameron & Truraiil d 197 lb . Combined population and physical 
evolution models were calculated for a representative sample of 
galaxies bv lEinastol ( 1972 , 1974 ). The last calculations showed 
that it was impossible to reproduce the rotation data by known 
stellar populations only. The only way to eliminate this conflict 
was to assume the presence of an unknown population - corona 
- with a very high value of the mass-to-light ratio, and a large 
radius and mass. Thus, the detailed modeling confirmed earlier 
results obtained by simpler models. But here we have one seri- 
ous difficulty - no known stellar population has so large a M/L 
value. 

Additional arguments for the presence of a spherical mas- 
sive population in spiral galaxies ca me from the stability criteri a 
against bar formation, suggested bv lOstriker & Peebles! dl973l) . 
Their numerical calculations demonstrated that initially very flat 
systems become rapidly thicker (during one revolution of the 
system) and evolve to a bar-like body. In real spiral galaxies a 
thin population exists, and it has no bar-like form. In their con- 
cluding remarks the authors write: "Presumably even Sc and 
other relatively 'pure' spirals must have some means of remain- 
ing stable, and the possibility exists that those systems also have 
very large, low -luminosity halos. The picture developed here 
agrees very well with the fact, noted b y several authors (see, 
for example, \Rogstad & Shostak (1972)), that the mass-to-light 
ratio increases rapidly with distance from the center in these 
systems; the increase may be due to the growing dominance 
of the high mass-to-light halo over the low mass-to-light ra- 
tio disk. It also suggests that the total mass of such systems 
has been severely und erestimated. In particular, the finding of 
Roberts & Rot s\nl 97^1 that the rotation curves of several nearby 



ground (CMB) radiation in combination with data from type la 
supernovae in nearby and very distant galaxies yield information 
on the curvature of the Universe that depends on the amount of 
Dark Matter and Dark Energy. 

Now we shall discuss these data in more detail. 

3.1 Stellar motions 

The local mass density near the Sun can be derived from vertical 
oscillations of stars n ear the galactic plane, as was discussed be 



fore. M odern data by Kuiiken & Gilmore ( 1989): lGilmore et al 



spirals become flat at large distances from the nucleus may in- 
dicate the presence of very extended halos having masses that 
diverge rapidly [M(r) prop to r] with distance." 

3 Dark Matter in astronomical data 

Modern astronomical methods yield a variety of independent in- 
formation on the presence and distribution of dark matter. For 
our Galaxy, the basic data are the stellar motions perpendicular 
to the plane of the Galaxy (for the local dark matter), the mo- 
tions of star and gas streams and the rotation (for the global dark 
matter). Important additional data come from gravitational mi- 
crolensing by invisible stars or planets. In nearby dwarf galax- 
ies the basic information comes from stellar motions. In more 
distant and giant galaxies the basic information comes from the 
rotation curves and the X-ray emission of the hot gas surround- 
ing galaxies. In clusters and groups of galaxies the gravitation 
field can be determined from relative motions of galaxies, the 
X-ray emission of hot gas and gravitational lensing. Finally, 
measurements of fluctuations of the Cosmic Microwave Back- 



(1989) have confirmed the results by Kuzmin and his collabora- 
tors. Thus we come to the conclusion that there is no evidence 
for the presence of large amounts of dark matter in the disk of 
the Galaxy. If there is some invisible matter near the galactic 
plane, then its amount is small, of the order of 15 percent of the 
total mass density. The local dark matter is probably baryonic 
(low-mass stars or jupiters), since non-baryonic matter is dissi- 
pationless and cannot form a highly flattened population. Spher- 
ical distribution of the local dark matter (in quantities suggested 
by Oort and Bahcall) is excluded since in this case the total mass 
of the dark population would be very large and would influence 
also the rotational velocity of the Galaxy at the location of the 
Solar System. 

Another information of the distribution of mass in the outer 
part of the Galaxy comes from streams of stars and gas. One of 
the streams discovered near the Galaxy is the Magellanic Stream 
of gas which forms a huge strip and c onnects the Large Magel- 



lanic Cloud (LMC) with the Galaxy (Mat hewson et all [1974). 



Model calculations emphasize that this stream is due to an en- 
counter of the LMC with the Galaxy. Kinematical data for the 
stream are available and support the hypoth esis on the presenc e 



of a massive halo surrounding the Galaxy dEinasto et al.L ll976a). 



Recently, streams of stars have been discovered within the 
Galaxy as well as around our giant neighbour M3 1 . Presently 
there are still few data on the kinematics of these streams. 

Several measurements of the dark mass halo were also per- 
formed using the motion of the satellite galaxies or the globular 
clusters. Measurements indicate the mass of the dark halo of 
about 2 x 1O I2 M . 

However, significant progress is expected in the near future. 
The astronomical satellite GAIA (to fly in 2011) is expected to 
measure distances and photometric data for millions of stars in 
the Galaxy. When these data are available, more information on 
the gravitation field of the Galaxy can be found. 

The motion of individual stars or gaseous clouds can be also 
studied in nearby dwarf galaxies. Determination of the dark halo 
was performed for over a dozen of them. Some of the newly 
discovered dwarfs, coming from the Sloan Digital Sky Survey, 
are very under-luminous but equally massive as the previously 
known dwarf galaxies in the Milky Way vicinity, which makes 
them good candidates for extreme examples of dark matter dom- 
inated objects. Also the studies of the disruption rate of these 
galaxies due to the interaction with the Milky Way imposes lim- 
its to the amount of dark mass in these objects. The results indi- 
cate that the dark matter in these systems exceeds by a factor a 
few the mass of stars. 




10-' 

Radius (kpc) 

Figure 2: The mean internal mass M(R) as a function of the 
radius R from the main galaxy in 105 pairs of galaxies (dots). 
Dashed line shows the contribution of visible populations, dot- 
ted line the contribution of the da rk corona, solid line the total 
distribution dEinasto et al l l!974ah . 
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Figure 3: Ma s ses ( in units 10 12 M Q ) of local giant galaxies 
(lOstriker et all 1 19741) (reproduced by permission of the AAS 



and authors). 



3.2 Dynamics and morphology of companion 
galaxies 

The rotation data available in early 1970s allowed to determine 
the mass distribution in galaxies up to their visible edges. In 
order to find how large and massive galactic coronas or halos 
are, more distant test particles are needed. If halos are large 
enough, then in pairs of galaxies the companion galaxies are 
located inside the halo, and their relative velocities can be used 
instead of the galaxy rotation velocities to find the distribution of 
m ass around giant galax ies. This test was made in dependently 
by lEinasto et alJ Jl974al) and lOstriker et all dl974h . see Figs. El 
and [3] The paper by Ostriker et al. begins with the statement: 
"There are reasons, increasing in number and quality, to believe 
that the masses of ordinary galaxies may have been underesti- 
mated by a factor of 10 or more" . The closing statement of the 
Einasto et al. paper is: "The mass of galactic coronas exceeds 
the mass of populations of known stars by one order of magni- 
tude. According to new estimates the total mass density of matter 
in galaxies is 20% of the critical cosmological density." 

The bottom line in both papers was: since the data suggest 
that all giant galaxies have massive halos/coronae, dark matter 
must be the dynamically dominating population in the whole 
Universe. 



Results of these papers were questioned bv lBurbidgd (1 1 9751) 
who noticed that satellites may be optical. To clari fy if the com 



panions are true members of the satellite systems, lEinasto et al 



(1974b) studied the morphology of companions. They found 
that companion galaxies are segregated morphologically: ellip- 
tical (non-gaseous) companions lie close to the primary galaxy 
whereas spiral and irregular (gaseous) companions of the same 
luminosity have larger distances from the primary galaxy; the 
distance of the segregation line from the primary galaxy de- 
pends on the luminosity of the satellite galaxy, Fig.|U This result 
shows, first of all, that the companions are real members of these 
systems - random by-fliers cannot have such properties. Second, 
this result demonstrates that diffuse matter has an important role 



in the evolution of galaxy systems. Morphological properties of 
companion galaxies can be explained, if we assume that (at least 
part of) the corona is gaseous. 

Additional arguments in favour of physical connection of 
companions with their primary galaxies came from the dynamics 
of small groups. Their mass distribution depends on the mor- 
phology: in systems with a bright primary galaxy the density 
(found from kinematical data) is systematically higher, and in el- 
liptical galaxy dominated systems it is also higher. The mass dis- 
tribution found from the kinematics of group members smoothly 
continues the mass distribution of the p rimary galaxies, found 



from rotation data ( lEinasto et al. , 1976b). 



3.3 Extended rotation curves of galaxies 

The dark matter problem was d iscussed in 1975 at two con 



ferences, in January in Tallinn dDoroshkevich et all 1 19751) and 
in July in Tbilisi. The central problems discussed in Tallinn 
were: Deuterium abun dance and the mean density of the uni- 



verse dZeldovichi 1 1 9751) . What is the physical nature of the dark 
matter? and: What is its role in the evolution of the Universe? 
Two basic models were suggested for coronas: faint stars or 
hot gas. It was found that both models have seriou s difficulties 
dJaaniste & SaarMl975l:lKomberg & Novikovlll975h . 

In Tbilisi the Third European Astronomical Meeting took 
place. Here the principal discussion was between the support- 
ers of the classical paradigm with conventional mass estimates 
of galaxies and of the new one with dark matter. The major ar- 
guments supporting the class ical paradigm were summarised by 
Mate rne & Tammannl d 1 976b . Their most serious argument was: 



Big Bang nucleosynthesis suggests a low-density Universe with 
the density parameter £2 « 0.05; the smoothness of the Hubble 
flow also favours a low-density Universe. 

It was clear that by sole discussion the presence and nature 
of dark matter cannot be solved, new data and more detailed 
studies were needed. The first very strong confirmation of the 
dark matter hypothesis came from new extended rotation curves 
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Figure 4: The distribution of luminosities of companion galaxies 
L(R) at various distances R from the main galaxy. Filled circles 
are for el liptical companions, o pen circles for spiral and irregular 
galaxies (lEinasto et al. . Il974bl) . 



of galaxies. 

In early 1970s optical data on rotation of galaxies were avail- 
able only for inner bright regions of galaxies. Radio observa- 
tions of the 21 -cm line reached much longer rotation curves well 
beyond the Holmberg ra dius of galaxies . All available rotation 
data were summarised by lRoberts ( 1975b in the IAU Symposium 
on Dynamics of Stellar Systems held in Besancon (France) in 
September 1974. Extended rotation curves were available for 14 
galaxies; for some galaxies data were available until the galacto- 
centric distance ~ 40 h Mpc (we use in this paper the Hubble 
constant in the units of Hq = 100 h km s Mpc -1 ), see Fig. Q] 
for M3 1 . About half of galaxies had flat rotation curves, the rest 
had rotation velocities that decreased slightly with distance. In 
all galaxies the local mass-to-light ratio in the periphery reached 
values over 100 in Solar units. To explain such high M/L values 
Roberts assumed that late-type dwarf stars dominate the periph- 
eral regions. 

In mid-1970s measurements of a number of spiral galaxies 
with the Westerbork Synthesis Radio Telescope were completed, 
and mass distr ibution mod els were built, all-together for 25 spi- 
ral galaxies dBosmalll978l) . see Fig. [5] Observations confirmed 
the general trend that the mean rotation curves remain flat over 
the whole observed range of distances from the center, up to 
~ 40 kpc for several galaxies. The internal mass within the ra- 
dius R increases over the whole distance interval. 

At the same time Vera Rubin and her collaborators devel- 
oped new sensitive detectors to measure optically the rotation 
curves of galaxies at very large galactocentric distances. Their 
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Figure 6: The integral masses as a function of the distance from 
the nucleus for sp iral galaxies of various morphological type 
( Rubin et al. , 1978) (reproduced by permission of the AAS and 
authors). 



results suggested that practi cally all spiral galaxies have ex- 
tended fiat rotation curves (IRubin et all 1 19781 Il980l) . The in- 
ternal mass of galaxies rises with distance almost linearly, up to 
the last measured point, see Fig. [6] 

These observational results confirmed the concept of the 
presence of dark halos of galaxies with a high confidence. 

Another very i mportant measurement was made by Faber 



and collaborators (Faber & Jackson! Il976t iFaber et all Il977t 



Faber & Gallagherl ll979h . They measured the central velocity 
dispersions for 25 elliptical galaxies and the rotation velocity of 
the Sombrero galaxy, a SO galaxy with a massive bulge and a 
very weak population of young stars and gas clouds just outside 
the main body of the bulge. Their data yielded for the bulge of 
the Sombrero galaxy a mass-to-light ratio M/L = 3, and for the 
mean mass-to-light ratios for elliptical galaxies about 7, close to 
the ratio for early type spiral galaxies. These observational data 
confirmed estimates based on the calculations of physical evolu- 
tion of galaxies, made under the assumption that the lower mass 
limit of the initial mass function (IMF) is for all galactic pop- 
ulations of the order of 0.1 M Q . These results showed that the 
mass-to-light ratios of stellar populations in spiral and elliptical 
galaxies are similar for a given colour, and the ratios are much 
lower than those accepted in earlier studies based on the dynam- 
ics of groups and clusters. In other words, high mass-to-light 
ratios of groups and clusters of galaxies cannot be explained by 
visible galactic populations. 

Earlier suggestions on the presence of mass discrepancy in 
galaxies and galaxy systems had been ignored by the astronom- 
ical community. This time new results were taken seriously. As 
noted by Kuhn, a scientific revolution begins when leading sci- 
entists in the field start to discuss the problem and arguments in 

favour of the new and the old paradigm. 

More data are slowly accumulating (ISofue & Rubinll200lh . 
New HI measurements from Westerbork extend the rotation 
curves up to 80 kpc (galaxy UGC 2487) or even 100 kpc 
(UGC 9133 and UGC 1 1852) showing flat rotation curves 



(INoordermeer et all 120051) . The HI distribution in the Milky 



Way has been r ecently studied up to distances of 40 kpc by 
Kalberlal d2003l) : iKalberla et all d2007l). The Milky Way rota- 
tion curve has been determined by IXue et al.1 (12008b up to ~ 60 



kpc from the study of ~ 2500 Blue Horizontal Branch stars from 
SDSS survey, and the rotation curves seems to be slightly falling 
from the 220 km s _1 value at the Sun location. Earlier determi- 
nations did not extend so far and extrapolations were affected 
by the presence of the ring-like structure in mass distribution at 
~ 14 kpc from the center. Implied values of the dark matter halo 
from different measurements still differ between themselves by 
a factor 2-3, being in the range from 10 12 - 2.5 x 1O 12 M . The 
central density of dark matter halos of galax i es is s urprisingly 
constant, about 0.1 M pc~ 3 dGilmore et all 12007b . Smallest 
dwarf galaxies have half-light radius about 120 pc, largest star 
clusters of similar absolute magnitude have half-light radius up 
to 35 pc; this gap sep arates systems with and without dark halos 



(Gilm oreetall 2008). 



3.4 X-ray data on galaxies and clusters of galax- 
ies 

Hot intra-cluster gas emitting X-rays was detected in almost all 
nearby clusters and in many groups of galaxies by the Einstein 
X-ray orbiting observatory. Observations confirmed that the hot 
gas is in hydrodynamical equilibrium, i.e. gas particles move in 
the general gravitation field of the clust er with velocities which 
correspond to the mass of the cl uster dForman & Jonesl U982; 



Sarazin, 1988; Rosati et al, 2002). 



The distribution of the mass in clusters can be determined 
if the density and the temperature of the intra-cluster gas are 
known. This method of determining the mass has a number of 
advantages over the use of the virial theorem. First, the gas is 
a collisional fluid, and particle velocities are isotropically dis- 
tributed, which is not true for galaxies as test particles of the 
cluster mass (uncertainties in the velocity anisotropy of galax- 
ies affect mass determinations). Second, the hydrostatic method 
gives the mass as a function of radius, rather than the total mass 
alone as given by the virial method. 

Using Einstein X-ray satellite data the method was applied 
to determine the mass of Coma, Perseu s and Virgo clusters 
(Bah call & Saraziri[l977tlMathewslll978h . The results were not 
very accurate since the temperature profile was known only ap- 
proximately. The results confirmed previous estimates of masses 
made with the virial method using galaxies as test particles. The 
mass of the hot gas itself is only about 0.1 of the total mass, 
approximately comparable to the luminous mass in galaxies. 

More recently clusters of galaxies have been observed in X- 
rays using the ROSAT satellite (operated in 1990-1999), and 
the XMM-Newton and Chandra observatories, launched both in 
1999. The ROSAT satellite was used to compile an all-sky cata- 
logue of X-ray clusters and galaxies. More than 1000 clusters up 
to a redshift ~ 0.5 were catalogued. Dark matter profi l es hav e 
been determined in a number of cases dHumphrev et al. , 2006). 

The XMM and Chandra observatories allow to get detailed 
images of X-ray cl usters, and to derive the density and te mpera- 
ture of the hot gas d Jordan et al. , 2004 ; Rasia et all 20061) . Using 
the XMM observatory, a survey of X-ray clusters was initiated to 
find a representative sample of clusters at redshifts up to z = 1 . 
The comparison of cluster properties at different redshifts allows 
to get more accurate information on the evolution of clusters 
which depend critically on the parameters of the cosmological 
model. 

Chandra observations allow to find the hot gas and total 
masses n ot only for group s and c lusters, but also for nearby 
galaxies (Hu mphrey et all (12006b . see also IMathews et al 



d2006l) . iLehmer et all d2008l) ). For early-type (elliptical) galax- 
ies the virial masses found were 0.7-9xl0 13 M . Local mass- 
to-light ratio profiles are flat within an optical half-light radius 
(R e ff), rising more than an order of magnitude at ~ lQR e ff, 
which confirms the presence of dark matter. The baryon frac- 
tion (most baryons are in the hot X-ray emitting gas) in these 
galaxies is fi, ~ 0.04 - 0.09. The gas mass profiles are simi- 
lar to the profiles of dark matter shifted to lower densities. The 
stellar mass-to-light ratios in these old bulge dominated galaxies 
are M*/L K ~ 0.5 -1.9 using the Salpeter IMF (for the infrared 
K-band, the ratios for the B-band are approximately 4 times 
higher). Interesting upper limits for the amount of hot plasma in 
the halo of the Milky Way were obtained in 2008 from the com- 
parative study of the tiny absorption lines in a few Galactic and 
extragalactic X-ray sources, giving the total column density of O 
VII less than 5 x 10 15 cirT 2 . Assuming that the gaseous baryonic 
corona has the mass of order of ~ 6 x 10 I0 M o , this measurement 
implies a very low metali city of the coron a plasma, below 3.7 
percent of the solar value dYao et al l l2008h . 



3.5 Galactic and extragalactic gravitational lens- 
ing 

Clusters, galaxies and even stars are so massive that their gravity 
bends and focuses the light from distant galaxies, quasars and 
stars that lie far behind. There are three classes of gravitational 
lensing: 

• Strong lensing, where there are easily visible distortions 
such as the formation of Einstein rings, arcs, and multiple 
images, see Fig. [7] 

• Weak lensing, where the distortions of background objects 
are much smaller and can only be detected by analyzing 
the shape distortions of a large number of objects. 

• Microlensing, where no shape distortion can be seen, but 
the amount of light received from a background object 
changes in time. The background source and the lens may 
be stars in the Milky Way or in nearby galaxies (M31, 
Magellanic Clouds). 







Figure 7: The Hubble Space Telescope image of the cluster 
Abell 2218. This cluster is so massive that its gravity bends 
the light of more distant background galaxies. Images of back- 
ground galaxies are distorted into stretched arcs (Astr. Pict. of 
the Day Jan. 11, 1998, Credit: W. Couch, R. Ellis). 

The strong lensing effect is observed in rich clusters, and 
allows to determine the distribution of the gravitating mass in 
clusters. Massive galaxies can distort images of distant single 



objects, such as quasars: as a result we observe multiple im- 
ages of the same quasar. The masses of clusters of galaxies de- 
termined using this method confirm the results obtained by the 
virial theorem and the X-ray data. 

Weak lensing allows to determine the distribution of dark 
matter in clusters as well as in superclusters. For the most lu- 
minous X-ray cluster known, RXJ 1347.5-1145 at the redshift 
z = 0.45, the lensing mass estimate is almost twice as high 
as that determined from the X-ray dat a. The mass-to-light ra- 
tio is MIL R = 200 + 50 in Solar units dFischer & Tvsoni ll997 



Fischer et all 11997). For other recent w o rk on weak len s ing an d 
X-ray clusters see iBradac et al.1 d2005l): [Dietrich et al.1 J2005); 



Clowe et all d2006bl) : Massey et all d2007h . 

A fraction of the invisible baryonic matter can lie in small 
compact objects - brown dwarf stars or jupiter-like objects. To 
find the fraction of these objects in the cosmic balance of matter, 
special studies have been initiated, based on the microlensing 
effect. 

Microlensing effects were used to find Massive Compact 
Halo Objects (MACHOs). MACHOs are small baryonic objects 
as planets, dead stars (white dwarfs) or brown dwarfs, which 
emit so little radiation that they are invisible most of the time. A 
MACHO may be detected when it passes in front of a star and 
the MACHOs gravity bends the light, causing the star to appear 
brighter. Several groups have used this method to search for the 
baryonic dark matter. Some authors claimed that up to 20 % of 
dark matter in our Galaxy can be in low-mass stars (white or 
brown dwarfs). However, observations using the Hubble Space 
Telescope's NICMOS instrument show that only about 6% of 
the stellar mass is composed of brown dwarfs. This corresponds 
to a negligible fraction of the total matter c ontent of the Universe 



dGraff & Freeselll99dlNa]ita et aUbOOOh 



4 The nature of Dark Matter 

By the end of 1970s most objections against the dark matter hy- 
pothesis were rejected. In particular, luminous populations of 
galaxies have found to have lower mass-to-light ratios than ex- 
pected previously, thus the presence of extra dark matter both 
in galaxies and clusters has been confirmed. However, there re- 
mained three problems: 

• It was not clear how to explain the Big Bang nucleosyn- 
thesis constraint on the low density of matter, and the 
smoothness of the Hubble flow. 

• If the massive halo (corona) is not stellar nor gaseous, of 
what stuff is it made of? 

• And a more general question: in Nature everything has its 
purpose. If 90 % of matter is dark, then there must be a 
reason for its presence. What is the role of dark matter in 
the history of the Universe? 

First we shall discuss baryons as dark matter candidates. 

4.1 Nucleosynthesis constraints on the amount of 
baryonic matter 

According to the Big Bang model, the Universe began in an ex- 
tremely hot and dense state. For the first second it was so hot 
that atomic nuclei could not form - space was filled with a hot 




baryonic matter is in the intergalactic matter, it is concentrated 
also in hot X-ray coronas of galaxies and clusters. 



Figure 8: The big-bang production of the light elements. The 
abundance of chemical elements is given as a function of 
the density of baryons, expressed in units of Q.bh 2 (horizon- 
tal axis). Predicted abundances are in agreement with mea- 
sured primeval abundances in a narrow range of baryon density 
dSchramm & Turnerlfl998l) . 



soup of protons, neutrons, electrons, photons and other short- 
lived particles. Occasionally a proton and a neutron collided 
and sticked together to form a nucleus of deuterium (a heavy iso- 
tope of hydrogen), but at such high tempera tures they were bro- 
ken i mmediately by high-energy photons ([Schramm & Turner , 
19981) . 

When the Universe cooled off, these high-energy photons 
became rare enough that it became possible for deuterium to 
survive. These deuterium nuclei could keep sticking to more 
protons and neutrons, forming nuclei of helium-3, helium-4, 
lithium, and beryllium. This process of element-formation is 
called "nucleosynthesis". The denser proton and neutron "gas" 
is at this time, the more of these light elements will be formed. 
As the Universe expands, however, the density of protons and 
neutrons decreases and the process slows down. Neutrons are 
unstable (with a lifetime of about 15 minutes) unless they are 
bound up inside a nucleus. After a few minutes, therefore, the 
free neutrons will be gone and nucleosynthesis will stop. There 
is only a small window of time in which nucleosynthesis can 
take place, and the relationship between the expansion rate of 
the Universe (related to the total matter density) and the density 
of protons and neutrons (the baryonic matter density) determines 
how much of each of these light elements are formed in the early 
Universe. 

According to nucleosynthesis data baryonic matter makes up 
0.04 of the critical cosmological density (Fig. [8). Only a small 
fraction, less than 10%, of the baryonic matter is condensed to 
visible stars, planets and other compact objects. Most of the 
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Figure 9: The evolution of mass-to-light ratios fs of ga lactic 
populations of different metal abundance Z ( Einastol 1972 ). The 
age of population t is expressed in years. An identical lower 
limit of IMF 0. 1 M was accepted. 



4.2 Baryonic Dark Matter 

Models of the galaxy evolution are based on stellar evolution 
tracks, star formation rates (as a f unction of time) , and the initial 
mass function (IMF). For IMF the Salpeter ( 19551) law is usually 
used: 

F(m) = am'", (1) 

where m is the mass of the forming star, and a and n are parame- 
ters. This law cannot be used for stars of arbitrary mass, because 
in this case the total mass of forming stars may be infinite. Thus 
we assume that this law is valid in the mass interval mo to m„ 
(the lower and upper limit of the forming stars, respectively). 

Early models o f physi cal e volution of gal axies were con- 
structed by Tinslev d 1968b and lEinasto ( 1972 ). These models 
show that the mass-to-light ratio M,/L, of the population i de- 
pends critically on the lower limit of the IMF, mo. It is natural to 
expect that in similar physical conditions (the metalicity of the 
gas in star formation regions) the lower mass limit of forming 
stars has similar values (Fig. [9]). An independent check of the 
correctness of the lower limit is provided by homogeneous stel- 
lar populations, such as star clusters. Here we can assume that 
all stars were formed simultaneously, the age of the cluster can 
be estimated from the HR diagram, and the mass derived from 
the kinematics of stars in the cluster. Such data are available for 
old metal-poor globular clusters, for relatively young medium- 
metal-rich open clusters as well as for metal-rich cores of galax- 
ies. This check suggests that in the first approximation for all 
populations similar lower mass limits (mo = 0.05 . . .0.1 M ) 
can be used; in contracting gas clouds above this limit the hy- 
drogen starts burning, below not. Using this mass lower limits 
we get for old metal -poor halo populations M, /L,- » 1, and for 
extremely metal-rich populations in central regions of galaxies 



Mil Lj = 10 . . . 30, as suggested by the central velocity disper- 
sion in luminous elliptical galaxies. For intermediate popula- 
tions (bulges and disks) one gets M,/L, = 3 ... 10, see Fig. [9] 
Modern data yield slightly lower values, due to more accurate 
measurements of velocity dispersions in the central regions of 
galaxi e s, as suggested in p ioneering studies by lFaber & JacksonI 
dl976h : lFaberetalJdl977L and more accurate input data for evo- 
lution models. 

To get very high values of M/L, as suggested by the dynam- 
ics of companion galaxies or rotation data in the periphery of 
galaxies, one needs to use a very small value of the mass lower 
limit ;«o <s 10~ 3 M . All known stellar populations have much 
lover mass-to-light values, and form continuous sequences in 
color- M/L and velocity dispersion-M/L diagrams. 

For this reason it is very difficult to explain the physical and 
kinematical properties of a stellar dark halo. Dark halo stars 
form an extended population around galaxies, and must have 
a much higher velocity dispersion than the stars belonging to 
the ordinary halo. No fast-movin g stars as possible can didates 
for stellar dark halos were found dJaaniste & Saaru 19751) . If the 
hypothetical population is of stellar origin, it must be formed 
much earlier than all known populations, because known stellar 
populations of different age and metalicity form a continuous 
sequence of kinematical and physical properties, and there is no 
place where to include this new population into this sequence. 
And, finally, it is known that star formation is not an efficient 
process - usually in a contracting gas cloud only about 1 % of 
the mass is converted to stars. Thus we have a problem how to 
convert, in an early stage of the evolution of the Universe, a large 
fraction of the primordial gas into this population of dark stars. 
Numerical simulations suggest, that in th early universe only a 
very small fraction of gas condenses to stars which ionize the 
rema i ning gas and stop for a ce rtain period further star formation 
(Eepil2003[lGao et allEo05bb . 

Stellar origin of dark matter in clusters was discussed by 
Napier & Guthrid (119751) ; they find that this is possible if the 
initial mass function of stars is strong l y bias ed toward very low- 
mass stars. iThorstensen & Partridge! d 19751) discussed the sug- 
gestion made by iTruran & Cameronl d 1 97 II) that there may have 
been a pre-galactic generation of stars (called now population 
III), all of them more massive than the Sun, which are now 
present as collapsed objects. They conclude that the total mass 
of this population is negligible, thus collapsed stars cannot make 
up the dark matter. 

Recently weak stellar halos have been detected around sev- 
eral nearby spiral galaxies at very large galactocentric distances. 
For instance, a ve ry weak stellar halo is found in M3 1 up t o dis- 
tance of 165 kpc dGilbert et all 120061: iKahrai et all 120061) . The 
stars of this halo have very low metalicity, but have anomalously 
red colour. The total luminosity and mass of these extended ha- 
los is, however, very small, thus these halos cannot be identified 
with the dark halo. 

Gaseous coron a s of g a laxie s and c l usters were discusse d 
in 1970s bv iFieldl dl972h. Isiikl dl974l) . iTarter & Silkl d 19741) . 
Romberg & NovikovT i 1975 ) and others. The general conclu- 
sion from these studies was that gaseous coronas of galaxies 
and clusters cannot consist of neutral gas since the intergalac- 
tic hot gas would ionise the coronal gas. On the other hand, a 
corona consisting of hot ionised gas would be observable. Mod- 
ern data show that part of the coronal matter around galaxies and 
in groups and clusters of galaxies consists indeed of the X-ray 
emitting hot gas, but the amount of this gas is not sufficient to 



explain the fiat rotation curves of galaxies (lTurneril2003l) . 

The result of these early discussions of the nature of dark 
halos were inconclusive - no appropriate candidate was found. 
For many astronomers this was an argument against the presence 
of dark halos. 

4.3 Non-baryonic Dark Matter and fluctuations 
of the CMB radiation 

Already in 1970s suggestions were made that some sort of non- 
baryonic elementary particles, such as massive neutrinos, mag- 
netic monopoles, axions, photinos, etc., may serve as candidates 
for dark matter particles. There were several reasons to search 
for non-baryonic particles as a dark matter candidate. First of 
all, no baryonic matter candidate did fit the observational data. 
Second, the total amount of dark matter is of the order of 0.2- 
0.3 in units of the critical cosmological density, while the nu- 
cleosynthesis constraints suggest that the amount of baryonic 
matter cannot be higher than about 0.04 of the critical density. 

A third very important observation was made which caused 
doubts to the baryonic matter as the dark matter candidate. In 
1964 Cosmic Microwave Background (CMB) radiation was de- 
tected. This discovery was a powerful confirmation of the Big 
Bang theory. Initially the Universe was very hot and all den- 
sity and temperature fluctuations of the primordial soup were 
damped by very intense radiation; the gas was ionized. But as 
the Universe expanded, the gas cooled and at a certain epoch 
called recombination the gas became neutral. From this time on 
density fluctuations in the gas had a chance to grow by gravi- 
tational instability. Matter is attracted to the regions were the 
density is higher and it flows away from low-density regions. 
But gravitational clustering is a very slow process. Model cal- 
culations show that in order to have time to build up all observed 
structures as galaxies, clusters, and superclusters, the amplitude 
of initial density fluctuations at the epoch of recombination must 
be of the order of 10 -3 of the density itself. These calculations 
also showed that density fluctuations are of the same order as 
temperature fluctuations. Thus astronomers started to search 
for temperature fluctuations of the CMB radiation. None were 
found. As the accuracy of measurement increased, lower and 
lower upper limits for the amplitude of CMB fluctuations were 
obtained. In late 1970s it was clear that the upper limits are 
much lower than the theoretically predicted limit 10~ 3 . 

Then astronomers recalled the possible existence of 
non-baryonic particles, such as heavy neutrinos. This sug- 
gestion was made in de penden tly by several as tro nomer s 
dCowsik & McCleUandl 
Tremaine & Gu nnl d 19791 
Chernii 



1973); Szalav & Mar: 



(1976); 
(1980b); 



»«. w vju.u, Doroshk evich et al.1 

nl dl98ll) : bond etall il983|)) and others. They found 



that, if dark matter consists of heavy neutrinos, then this helps 
to explain the paradox of small temperature fluctuations of the 
cosmic microwave background radiation. This problem was 
discussed in a conference in Tallinn in April 1981. Recent 
experiments by a Moscow physicist Lyubimov were announced, 
which suggested that neutrinos have masses. If so, then the 
growth of perturbations in a neutrino-dominated medium can 
start much earlier than in a baryonic medium, and at the time of 
recombination perturbations may have amplitudes large enough 
for structure formation. The Lyubimov results were never con- 
firmed, but it gave cosmologists an impulse to take non-baryonic 
dark matter seriously. In the conference banquet Zeldovich 
gave an enthusiastic speech: "Observers work hard in sleepless 



nights to collect data; theorists interpret observations, are often 
in error, correct their errors and try again; and there are only 
very rare moments of clarification. Today it is one of such rare 
moments when we have a holy feeling of understanding the 
secrets of Nature." Non-baryonic dark matter is needed to start 
structure formation early enough. This example illustrates well 
the attitude of theorists to new observational discoveries - the 
Eddington's test: 'Wo experimental resu lt should be b elieved 
until confirmed by theory " (cited by Mike lTurnerl d2000l) ). Dark 
matter condenses at early epoch and forms potential wells, 
the baryonic matter fl ows into these wells and forms galaxies 
dWhite&Reeslll978l) . 

The search of dark matter can also be illustrated with the 
words of Sherlock Holmes "When you have eliminated the im- 
possible, whate ver remains, however impr obable, must be the 
truth" (cited bv lBinnev & Tremaind ( ll987l) ). The non-baryonic 
nature of dark matter explains the role of dark matter in the evo- 
lution of the Universe, and the discrepancy between the total 
cosmological density of matter and the density of baryonic mat- 
ter, as found from the nucleosynthesis constraint. Later studies 
have demonstrated that neutrinos are not the best candidates for 
the non-baryonic dark matter, see below. 

4.4 Alternatives to Dark Matter 

The presence of large amounts of matter of unknown origin 
has given rise to speculations on the validity of the New- 
ton law of gravity at very large distances. One of such at- 
temp ts is the Modified Newton D ynamics (MOND), sug gested 
by M ilgrom & Bekensteinl (119871) . for a discussion see ISanders 
(1990)." Indeed, MOND is able to represent a number of ob- 
servational data without assuming the presence of some hidden 
matter. However, there exist several arguments which make this 
model unrealistic. 

First of all, in the absence of large amounts of non-baryonic 
matter during the radiation domination era of the evolution of the 
Universe it would be impossible to get for the relative amplitude 
of density fluctuations a value of the order of 1CT 3 , needed to 
form all observed structures. 

Second, there exist numerous direct observations of the dis- 
tribution of mass, visible galaxies and the hot X-ray gas, which 
cannot be explained in the MOND framewo rk. One of such ex 



amples is the "bullet" cluster IE 0657-558 dClowe et al.L 12004 



Markevitch et all 12004 Iciowe et all l2006al) . shown in Fig. |T0j 
This is a pair of galaxy clusters, where the smaller cluster (bul- 
let) has passed the primary cluster almost tangentially to the line 
of sight. The hot X-ray gas has been separated by ram pressure- 
stripping during the passage. Weak gravitation lensing yields 
the distribution of mass in the cluster pair. Lensing observations 
show that the distribution of matter is identical with the distribu- 
tion of galaxies. The dominant population of the baryonic mass 
is in X-ray gas which is well separated from the distribution of 
mass. This separation is only possible if the mass is in the col- 
lisionless component, i.e. in the non-baryonic dark matter halo, 
not in the baryonic X-ray gas. 

5 Dark Matter and structure formation 

It is clear that if dark matter dominates in the matter budget of 
the Universe, then the properties of dark matter particles deter- 
mine the formation and evolution of the structure of the Uni- 



verse. In this way the dark matter problem is related to the large- 
scale structure of the Universe. 

5.1 The distribution of galaxies and clusters 

Already in the New General Catalogue (NGC) of nebulae, com- 
posed from observations by William and John Herschel, a rich 
collectio n of nearby galaxies in the Virgo constellation was 



known. Ide Vaucouleursl (I1953al) called this system the Local 



Super-galaxy, presently it is known as the Virgo or Local Su- 
percluster. Detailed investigation of the distribution of galaxies 
became possible when Harlow Shapley started in the Harvard 
Observatory a systematic photogra phic survey of galaxies in se- 
lected areas, up to 18th magnitude (Shapley, 1940). Shapley dis- 
covered several other rich superclusters, one of them is presently 
named the Shapley Supercluster. These studies showed also that 
the mean spatial density of galaxies is approximately indepen- 
dent of the distance and of the direction in the sky. In other 
words, the Harvard survey indicated that galaxies are distributed 
in space more-or-less homogeneously, as expected from the gen- 
eral cosmological principle. 

A complete photographic survey of galaxies was made in 
the Lick Observatory wit h the 20-inch Carnegie astrograph by 
Shane & Wirtanenl ( 1967 ). Galaxy counts were made in cells 
of size 10' x 10', and the distribution of the number density of 
galaxies was studied. The general conclusion from this study 
was that galaxies are mostly located in clusters, the number of 
galaxies per cluster varying widely from pairs to very rich clus- 
ters of the Coma cluster type. The Lick counts were reduced by 
Jim Peebles and collaborators to exclude count limit irregulari- 
ties; the resulting distribution of galaxies in the sky is shown in 
Fig-ttU 

A much deeper photographic survey was made using the 48- 
inch Palomar Schmidt telescope. Fritz Zwicky used this survey 
to compile for the North ern hemisphere a ca talogue of galaxies 



and clusters of galaxies (Zw icky et allll968l) . The galaxy cata- 
logue is complete up to 15.5 photographic magnitude. George 
Abell used the same survey to compile a catalogue of rich clus- 
ters of galaxies for the North ern sky ? later the catalogue wa s 
continued to the Southern sky dAbellLll958l:lAbell et allf!^ ). 
Using apparent magnitudes of galaxies approximate distances 
(distance classes) were estimated for clusters in both catalogues. 
Authors noticed that clusters of galaxies also show a tendency of 
clustering, similar to galaxies which cluster to form groups and 
clusters. Abell called these second order clusters superclusters, 
Zwicky - clouds of galaxies. 

The Lick counts as well as galaxy and cluster cata- 
logues by Zwicky and Ab e ll wer e analysed by Jim Peeble s 
and collaborators jPeeblesI (fl973l) lHauser & Peebles! dl973h . 
Peebles & Hauseri dl974 . |PeeblesT(ll974l) ). To describe the dis- 
tribution of galaxies Peebles introduced the two-point corre- 
lation (or covariance) function of galaxies ^Peebles & GrothL 
1975UGroth & Peebles! 1 19771 iFry & PeeblesLll978l) . This func- 
tion describes the probability to find a neighbour at a given an- 
gular separation in the sky from a galaxy. At small separations 
the spatial galaxy correlation function can be approximated by 
a power law: ^ = (r/ro)~ r , with the index y = 1.77 + 0.04. 
The distance ro, at which the correlation function equals unity, 
is called the correlation length. For galaxy samples its value is 
ro ~ 5 hT l Mpc, and for clusters of galaxies ro ~ 30 /T 1 Mpc. On 
scales > 5 times the correlation length the correlation function 
is very close to zero, i.e. the distribution of galaxies (clusters) is 
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Figure 10: Images of the merging 'bullet' cluster 1E0657-558. The left panel shows a direct image of the cluster obtained with 
the 6.5-m Magellan telescope in the Las Campanas Observatory, the right panel is a X-ray satellite Chandra image of the cluster. 
Shock waves of the gas are visible, the gas of the smaller 'bullet' cluster (right) lags behind the cluster galaxies. In both panels 
green contours are equidensity levels of the gravitational potential of the cluster, found using weak gravitational lensing of distant 
galaxies. The white bar has 200 kpc/h length at the distance of the cluster. Note that contours of the gravitational potential coincide 
with the location of visible galaxies, but not with the location of the X-ray gas (the dominant baryonic component of clusters) 
( Clowe et al. . 2006a) (reproduced by permission of the AAS and authors). 



essentially random. 

The conclusion from these studies, based on the appar- 
ent (2-dimensional) distribution of galaxi es and cluste rs in 
the sky confirmed the picture suggested by iKiangl dl967h and 
de Vaucouleursl ( 1970h . among others, that galaxies are hierar- 
chically clustered. However, this hierarchy does not continue to 
very large scales as this contradicts observations, which show 
that on very large scales the distribution is homogeneous. A 
theoretical explanation of this picture was given by Peebles 
in his hierarchical c lustering scenario of structure formation 
dPeebles & YuLll970UPeebleslll97ll) . 



5.2 Superclusters, filaments and voids 

In 1970s new sensitive detectors were developed which allowed 
to measure redshifts of galaxies up to fainter magnitudes. Tak- 
ing advance of this development several groups started to in- 
vestigate the environment of relatively rich clusters of galax- 
ies, such as the Coma cluster and clusters in the Hercules su- 
perc luster, with a limiting magnitude about 15.5. During this 
study Chincarini, Gregory, Rood, Thompson and Tifft noticed 
that the main clusters of the Coma supercluster, A 1656 and 
A 1367, are surrounded by numerous galaxies, forming a cloud 
around clusters at the redshift ~7000 km/s. The Coma super- 
cluster lies behind the Local supercluster, thus galaxies of the 
Local supercluster also form a condensation in the same direc- 
tion at the redshift about 1000 km/s. In between there is a group 
of galaxies around NGC 4169 at the redshift ~4000 km/s, and 
the s pace between these systems is completely devoid of ga lax- 
ies dChincarini & Roodlll976l : lGregorv & Thompsonill978l) . A 



similar picture was observed in front of the Hercules and Perseus 
superclusters. 

In 1970s there were two mai n rivaling theories of structure 
formation: the pancake theory by Zeldovich] (| 1 970» . and the hi- 
erarchical clustering theory by iPeeblesI d 197 lb . According to 
the Zeldovich scenario the structure forms top-down: first mat- 
ter collects into pancakes and then fragments to form smaller 
units. In the hierarchical clustering scenario the order of the 



formation of systems is the opposite: first small-scale systems 
(star-cluster sized objects) form, and by clustering systems of 
larger size (galaxies, clusters of galaxies) form; this is a bottom- 
up scenario. 

Zeldovich asked Tartu astronomers for help in solving the 
question: Can we find observational evidence which can be used 
to discriminate between various theories of galaxy formation? In 
solving the Zeldovich question we started from the observational 
fact suggesting that random velocities of galaxies are of the or- 
der of several hundred km/s. Thus during the whole lifetime 
of the Universe galaxies have moved from their place of origin 
only by about 1 h Mpc. In other words - if there exist some 
regularities in the distribution of galaxies, then these regularities 
must reflect the conditions in the Universe during the formation 
of galaxies. 

In mid-1970s first all-sky complete redshift surveys of g alax- 
ies were just available: the Ide Vaucouleurs et alj (11976b Sec- 
ond Revised Catalogue of Galaxies, the Sh apley- Adams revised 
catalogue by ISandage & Tammannl dl98ll) . complete up to the 
magnitude 13.5 (new redshifts were available earlier (Sandage, 
19781) ). The common practice to visualise the three-dimensional 
distribution of galaxies, groups and clusters of galaxies is the 
use of wedge-diagrams. In these diagrams, where galaxies as 
well as groups and clusters of galaxies were plotted, a regu- 
larity was clearly seen: galaxies and clusters are concentrated 
to identical essentially one-dimensional syst ems, and the space 
betwe en these systems is practically empty dJoeveer & Einasto, 
19781) . This distribution was quite similar to the distribution 



of test particles in a numerical simulation of the evolution of 
the structure of the Universe prepared by the Zeldovich group 
dDoroshkevich et al.l dl980al) . early results of simulation were 
available already in 1976). In this simulation a network of high- 
and low-density regions was seen: high-density regions form 
cells which surround large under-dense regions. Thus the ob- 
served high-density regions could be identified with Zeldovich 
pancakes. 

The Large-Scale Structure of the Universe was discussed 
at the IAU symposium in Tallinn 1977. The amazing prop- 



Figure 11: The two-dimensional distribution of galaxies according to the Lick counts (ISeldner et all 119771) . The north galactic 
pole is at the center, the galactic equator is at the edge. Superclusters are well seen, the Coma cluster is located near the center 
(reproduced by permission of the AAS and authors). 



erties of the distributio n of galaxies w ere reported by four 
differen t groups: iTuUv & Fished (119781) for the Local super- 
cluster, Joeveer & Einastd (119781) for the Perseus superclus- 
ter, Tarenghi et al. d 19781) for the Hercules supercluster, and 
Tifft & Gregory ( 19781) for th e Perseus supercluster; see also 
Gregory & Thompsonl ( 1978 ) for the Coma supercluster and 
Joeve er et al dl978h for the distribution of galaxies and clus- 
ters in the Southern galactic hemisphere. The presence of voids 
(holes) in galaxy distribution was suggested in all four reports. 
Tully & Fisher demonstrated a movie showing a filamentary 
distribution of galaxies in the Local supercluster. Joeveer and 
Einasto emphasized the presence of fine structure: groups and 
clusters of galaxies form chains in superclusters and connect su- 
perclusters to an infinite network, as seen from wedge diagrams 
in Fig. [12] They demonstrated also morphological properties 
of the structure of superclusters: clusters and groups within the 
chain are elongated along the chain, and main galaxies of clus- 
ters (supergiant galaxies of type cD) are also elongated along 
the chain. A long chain of clusters, groups and galaxies of the 
Perseus-Pisces supercluster is located almost perpendicular to 
the line of sight. The scatters of positions of clusters/groups 



along the chain in the radial (redshift) and tangential directions 
are practically identical. This demonstrates that the chain is es- 
sentially an one-dimensional structure. 

A direct consequence from this observation is that galaxies 
and groups/clusters of the chain are already formed within the 
chain. A later inflow from random locations to the chain is ex- 
cluded, since in this case it would be impossible to stop galaxies 
and clusters in the chain after the inflow. The main results of the 
symposium were summarised by Malcolm Longair as follows: 
To me, some of the most exiting results presented at this sym- 
posium concerned the structure of the Universe on the largest 
scales. Everyone seemed to agree about the existence of super- 
clusters on scales ~ 30 — 100 Mpc. But perhaps even more 
surprising are the great holes in the Universe. Peebles' pic- 
ture, Einasto 's analysis of the velocity distribution of galaxies 
which suggests a "cell-structure" and Tiffts's similar analysis 
argue that galaxies are found in interlocking chains over scales 
~ 50 - 100 Mpc forming pattern similar to a lace-tablecloth. 

New data gav e strong support to the pancake scenario by 
Zeldovichl dl978l) . However, some important differences be- 
tween the model and observations were evident. First of all, 
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Figure 12: Wedge diagrams for two declination zones. Filled circles show rich clusters of galaxies, open circles - groups, dots - 
galaxies, crosses - Markarian galaxies. In the 15° - 30° zone two rich clusters at RA about 12 h are the main clusters of the Coma 
supercluster, in the 30° - 45° zone clusters at RA about 3 h belong to the main chain of clusters and galaxies of the Perseus-Pisces 
supercluster. Note the complete absence of galaxies in front of the Pers eus-Pisces supercluster, an d galaxy chains leading from the 
Local supercluster towards the Coma and Perseus-Pisces superclusters (iJoeveer & EinastoUl978l) . 
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Figure 13: A sli ce of the Universe accordi ng to the CfA Second 
redshift survey (Ide Lapparent et all 1 1986b . Galaxy chains con- 
necting the Local and Coma superclusters are seen more clearly; 
the connection between the Hercules and Coma superclusters is 
also visible (reproduced by permission of the AAS and authors). 



numerical simulations showed that there exists a rarefied pop- 
ulation of test particles in voids absent in real data. This was 
the first indication for the presence of physical biasing in galaxy 
formation - there is primordial gas and dark matter in voids, but 
due to low density no galaxy formation takes place here. Theo- 
retical explanat ion of the ab s ence o f galaxies in voids was given 



by Enn Saar ( Einasto et al. . 1980h . In over-dense regions the 



density increases until the matter collapses to form compact ob- 
jects (Zeldovich pancakes). In under-dense regions the density 
decreases exponentially but never reaches a zero value - gravity 
cannot evacuate voids completely. 

The second difference lies in the structure of galaxy sys- 
tems in high-density regions: in the original pancake model 
large-scale structures (superclusters) have rather diffuse forms, 
real superclusters consi st of multiple intertwined filaments : 
Joeveer & Einastol dl978l) . IZeldovich etail dl982h . lOortl dl983l) . 
In the original pancake scenario small-scale perturbations were 



damped. This scenario corresponds to the neutrino-dominated 
dark matter. Neutrinos move with very high velocities which 
wash out small-scale fluctuations. Also, in the neutrino- 
dominated Universe superclusters and galaxies within them 
form relatively late, but the age of old stellar populations in 
galaxies suggests an early start of galaxy formation, soon af- 
ter the recombination epoch. In other words, the original 

n: — : — — n r — — 1 



panca ke scenario was in tr ouble (Bondetal., 1982; Peebles 
1982j: IZeldovich et all[l98llBond & SzalayLll983HWhite et al 



19831). 

The presence of voids in galaxy distribution was initially 
met with scepticism, since 3-dimensional data were available 
only for bright galaxies, and faint galaxies could fill voids. 
However, independent evidence w as soon found. A very large 
void was discovered in Bootes by iKirshner et al.l d 198 lb . The 
filamentary nature of galaxy distribution is very clearly seen 
in the 2nd Center for Astrophysics ( Harvard) Redshift Surve 



by Huchra, Geller a nd collaborators dde Lapparent et all 1198 



Huchra et all 119881) . complete up to 15.5 apparent blue magni- 



tude in the Northern Galactic hemisphere, see Fig. IT3l 

Huchra initiated a near-infrared survey of nearby galax- 
ies, the Two M i cron A ll-Sky Survey (2MASS) dHuchral l2000t 
2006). Photometry in 3 near-infrared spectral 



Skrutskie et al 



bands is completed, it includes about half a million galaxies up 
to the limiting K magnitude 13.5. The redshifts are planned to be 
measured for all galaxies up to K = 1 1 .25. The advantage of this 
survey is the coverage of low galactic latitudes up to 5 degrees 
from the Galactic equator. For the Southern sky the redshift sur- 
vey of 2MASS galaxies is almost completed using the 6 degree 
Field Survey with the Australian large Schmidt telescope. The 
filamentary character of the distribution of galaxies is very well 
seen. 

A much deeper redshift survey up to the blue magnitude 
19.4 was recently completed using the Anglo-Australian 4- 
m telescope. This Two degree Field Galaxy Redshift Survey 
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Figure 14: A wedge diagram of the luminosity density field of galaxies of the 2dFGRS Northern equatorial zone, +1.5 degrees 
around the equator. The luminosity densities have been corrected for the incompleteness effect; the RA coordinate is shifted 
so that the plot is symmetrical around the vertical axis. The rich superc luster at the distance ~ 250 hr x Mpc from the observer is 
SCL126, according to the supercluster catalogue by lEinasto et al.l (120011) . This figure illustrates the structure of the cosmic web (the 
supercluster-void network). The density field shows that rich superclusters contain many rich clusters of galaxies, seen in picture 



as red dots. Filaments consisting of less luminous galaxies and located between superclusters are also clearly seen (lEinasto et al 
20071). 



(2dFGRS) covers an equatorial strip in the Northern Galactic 
hemisphere and a contiguous area in t he Southern hemisphere 



(Coll ess et aUBOOlUCross et al.L l2001). Over 250 thousand red 



shifts have been measured, which allows to see and measure the 
cosmic web (supercluster-void network) up to the redshift 0.2, 
corresponding to a co-moving distance about 575 h' 1 Mpc. The 
luminosity density field calculated for the Northern equatorial 
slice of the 2dFGRS is shown in Fig. [14] 

Presently the largest project to map the Universe, the Sloan 
Digital Sky Survey (SDSS) mentioned already before, has been 
initiated by a number of Am erican, Japanese and European uni- 
versities and observatori es (lYork et al. , 2000; IStoughton et al 



2002; IZehavi et al.Ll2002h . The goal is to map a quarter of the en 
tire sky: to determine positions and photometric data in 5 spec- 
tral bands of galaxies and quasars of about 100 million objects 
down to the red magnitude r = 23, and redshifts of all galaxies 
down to r = 17.7 (about 1 million galaxies), as well as the red- 
shifts of Luminous Red Giants (LRG, mostly central galaxies of 



groups and clusters) down to the absolute magnitude about -20. 
All 7 data releases have been made public. This has allowed 
to map the largest volume of the Universe so far. LRGs have a 
spatial density about 10 times higher than rich Abell clusters of 
galaxies, which allows to sample the cosmic web with sufficient 
details up to a redshift ~ 0.5. 

5.3 Structure formation in the Cold Dark Matter 
scenario 



A consistent picture of the structure formation in the Universe 
slowly emerged from the advancement in the observational stud- 
ies of the large scale structure and the Cosmic Microwave Back- 
ground, and in the development of theory. The most important 
steps along the way were the following. 

Motivated by the observational problems with neu- 
trino dark matter, a new d ar k mat ter scenario w a s sug - 
gested by iBlumenthal etaD dl982h : iBond et al.1 dl982h : 



IPagels & Primackl dl982h: IPeebles dl982b : iBond & Szalav 
dl983b : boroshkevich & Khlopovl dl984l) with hypothetical 



particles as axions, gravitinos, photinos or unstable neutrinos 
playing the role of dark matter. This model was called the Cold 
Dark Matter (CDM) model, in contrast to the neutrino-based 
Hot Dark Matter model. Newly suggested dark matter particles 
move slowly, thus small-scale perturbations are not suppressed, 
which allows an early start of the structure formation and the 
formation o f fine structure. Advan tages of this model were 
discussed bv lBlumenthal et alj dl984 . 

Next, the cosmological constant, A, was incorporated into 
the scheme. Arguments favouring a m odel with the cosmologi - 
cal constant were s u ggested already by Gunn & Tinslevl d 19751) : 



Turn er et al.l (11984); Ko fman & Starobi nskii (1985): combined 



constraints on the density of the Universe, ages of galaxies, and 
baryon nucleosynthesis. 

Finally, there was a change in the understanding of the for- 
mation of initial perturbations which later lead to the observed 
structure. To explain the fl atness of the Universe the inflation 
scenario was suggested by iGuthl d 198 lb . iKofman et al.l d 1985b 
and others. According to the inflation model in the first stage the 
expansion of the Universe was extremely rapid (this early stage 
is called the inflation epoch). Such an evolutionary scenario al- 
lows for the creation of the visible part of the Universe out of a 
small causally connected region and explains why in the large 
scale the Universe seems roughly uniform. Perturbations of the 
field are generated by small quantum fluctuations. These pertur- 
bations form a Gaussian random field, they are scale-invariant 
and have a purely adiabatic primordial power spectrum. 

These ideas were progressively incorporated into the com- 
puter simulations of increasing complexity. 

Pioneering numerical simulations of the evolution of 
the structure of the Universe were made in 1970s by 



Milled dl978h. lAarseth et al.l d 1979b and the Zeldovich group 
dDoroshkevich et al.Lll980al) 7 using direct numerical integration. 
In early 1980s the Fourier transform was suggested to calculate 
the force field which allowed to increase the number of test par- 
ticles. 

Numerical simulations of structure evolution for the hot and 



cold dark m a tter were co mpared by lMelott et al.l dl983l) . and by 



White et alj d!983l 11987b (standard CDM model with density 
parameter Q,„ = 1). In contrast to the HDM model, in the CDM 
scenario the structure formation starts at an early epoch, and su- 
perclusters consist of a network of small galaxy filaments, simi- 
lar to the observed distribution of galaxies. Thus CDM simula- 
tions reproduce quite well the observed structure with clusters, 
filaments and voids, including quantitative characteristics (per- 
colation or connectivity, the multiplicity distribution of systems 
of galaxies). 

Model s with the cosmological A-term were developed by 



Gramann ( 1988). Comparison of the SCDM and ACDM models 



shows that the structure of the cosmic web is similar in both 
models. However, in order to get a correct amplitude of density 
fluctuations, the evolution of the SCDM model has to be stopped 
at an earlier epoch. 

The largest so far simulation of the evolution of the 
structure - the Millennium Simulation - was made in the 
Max-Planck Institute for Ast rophysics in Garching by Volker 



Spring el and collaborators dSpringel et all 120051 : iGao et al 



l2005d: ISpringel et all EooHb . The simulation is assuming the 



particles that allowed to follow the evolution of small-scale fea- 
tures in galaxies. Using a semi-analytic model the formation and 
evolution of galax i es was also simulate d dDi Matteo et all [2005; 
Gao et al. , 2005b; ICroton et alj , 2006). For simulated galaxies 
photometric properties, masses, luminosities and sizes of princi- 
pal components (bulge, disk) were found. The comparison of 
this simulated galaxy catalogue with observations shows that 
the simulation was very successful. The results of the Millen- 
nium Simulation are frequently used as a starting point for fur- 
ther more detailed simulations of evolution of single galaxies. 

One difficulty of the original pancake scenario by Zeldovich 
is the shape of objects formed during the collapse. It was 
assumed that forming systems are flat pancake-like objects, 
whereas dominant features of the cosmic we b are filamen t s. Thi s 
discrepancy has found an explanation by IBond et al.l 0996). 
They showed that in most cases just essentially one-dimensional 
structures, i.e. filaments form. 

The ACDM model of structure formation and evolution 
combines essential aspects of both original structure formation 
models, the pancake and the hierarchical clustering scenario. 
First structures form at very early epochs soon after the recom- 
bination in places where the primordial matter has the highest 
density. This occurs in the central regions of future superclus- 
ters. First objects to form are small dwarf galaxies, which grow 
by infall of primordial matter and other small galaxies. Thus, 
soon after the formation of the central galaxy other galaxies fall 
into the gravitational potential well of the supercluster. These 
clusters have had many merger events and have "eaten" all its 
nearby companions. During each merger event the cluster suf- 
fers a slight shift of its position. As merger galaxies come from 
all directions, the cluster sets more and more accurately to the 
center of the gravitational well of the supercluster. This explains 
the fact that very rich clusters have almost no residual motion in 
respect to the smooth Hubble flow. Numerous examples of the 
galaxy mergers are seen in the images of galaxies collected by 
the Hubble Space Telescope. 



6 Matter-energy content of the Universe 
6.1 Dark Matter and Dark Energy 

In early papers on dark matter the total density due to visible 
and dark matter was estimated to be about 0.2 of the critical 
cosmological density. These estimates were based on the dy- 
namics of galaxies in groups and clusters. This density estimate 
can be interpreted in two different ways: either we live in an 
open Universe where the total density is less than the critical 
density, or there exists some additional form of matter/energy 
which allows the Universe to be closed, i.e. to have the critical 
total density. The additional term was identified with the Ein- 
stein A-term, so that the total matter/energy de nsity was taken to 
be equal to the critical cosmological density dG unn & Ti nslev , 



1975; Turner et al., 1984; Kofman & Starobinskii, 1985). Ini 



ACDM initial power spectrum. A cube of the comoving size of 
500 hT x Mpc was simulated using about 10 billion dark matter 



tially there was no direct observational evidence in favour of 
this solution and it was supported basically on general theoret- 
ical grounds. In its early evolution the size of the Universe in- 
creases very rapidly and any deviation from the exact critical 
density would lead to a rapid change of the relative density, ei- 
ther to zero, if the initial density was a bit less than the critical 
one, or to infinity, if it was greater than critical. In other words, 
some fine tuning is needed to keep the density at all times equal 
to the critical one. 
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Figure 15: Upper panel shows the acoustic peaks in the angular 
power spectrum of the CMB radiation according to the WMAP 
and other recent data, compared with the ACDM model using all 
available data. The lower panel shows the signature of baryonic 
acoustic oscillations in t he matter t wo-point correlation function 
( Eisenstein et al. , 2005 : Kolbl 2007 ) (reproduced by permission 
of the author) 



In subsequent years several new independent methods were 
applied to estimate the cosmological parameters. Of these new 
methods two desire special attention. One of them is based 
on the measurements of small fluctuations of the Cosmic Mi- 
crowave Background (CMB) radiation, and the other on the ob- 
servation of distant supernovae. 

According to the present cosmological paradigm the Uni- 
verse was initially very hot and ionized. The photons provided 
high pressure and prevented baryons from moving. Perturba- 
tions of baryons did not grow, but oscillated as sound waves. 
The largest possible wavelength of these oscillations is given 
by the sound horizon size at the decoupling. This wavelength 
is seen as the first maximum in the angular power spectrum 
of the CMB radiation. The following maxima correspond to 
overtones of the first one. The fluctuations of CMB radiation 
were first detected by the COBE satellite. The first CMB data 
were not very accurate, since fluctuations are very small, of the 
order of 10~ 5 . Subsequent experiments carried out using bal- 
loons, ground based instruments, and more recently the Wilkin- 
son Microwave Anisotropy Probe (WMAP) satellite, allowed 
to measure the CMB ra diation and it s powe r spectrum with a 



much higher precision ( Spergel et all 120031) . The position of 



matter/energy density. Observations confirm the theoretically 
favoured value 1 in units of the critical cosmological density, 
see Fig.[T31 

When recombination begins, the small overdensities of 
baryon gas launch spherical shock waves in the photon-baryon 
fluid. After some time photons completely decouple from 
baryons, and the baryons loose photon pressure support. The 
shock stops after traveling a distance of about 150 Mpc (in co- 
moving coordinates). This leads to an overdensity of the bary- 
onic medium on a distance scale of 150 Mpc. This overdensity 
has been recently detected in the correlation f unction of Lumi- 
nous Red Giant ga laxies of the SDSS survey feisenstein et all 
2005;|Hutsl|2b06), see lower panel of Fig. [15] Baryonic acous- 
tic oscillations depend on both the total matter/energy density 
and the baryon density, thus allowing to estimate these parame- 
ters. 

Another independent source of information on cosmologi- 
cal parameters co mes from the distant sup ernova experiments. 
Two teams, led bv lRiess et alj dl998ll2007t) (High-Z Supernova 
Search Team) and Perlmutte r et al. Hi 9991) (Supernova Cosmol- 
ogy Project), initiated programs to detect distant type la super- 
novae in the early stage of their evolution, and to investigate 
with large telescopes their properties. These supernovae have 
an almost constant intrinsic brightness (depending slightly on 
their evolution). By comparing the luminosities and redshifts of 
nearby and distant supernovae it is possible to calculate how fast 
the Universe was expanding at different times. The supernova 
observations give strong support to the cosmological model with 
the A term, see Fig. [T6l 
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Figure 16: Results of the Supernova Legacy Survey: appar- 
ent magnitudes of supernovae are normalised to the standard 
ACDM model, shown as solid line. Dashed line shows the 



Einstein-de Sitter model with Q m = 1 dKolbi 12007 ) (reproduced 
by permission of the author). 

Different types of dark energy affect the rate at which the 
Universe expands, depending on their effective equation of state. 
The cosmological constant has one equation of state. The other 
possible candidate of dark energy is quintessence (a scalar field) 
that has a different equation of state. Each variant of dark energy 
has its own equation of state that produce s a signature in the 
Hubble diagram of the type la supernovae (ITurnen.120031) . 

The combination of the CMB and supernova data allows to 
estimate independently the matter density and the density due to 
dark energy, shown in Fig. [17] The results of this combined ap- 
proach imply that the Universe is expanding at an accelerating 
rate. The acceleration is due to the existence of some previ- 
ously unknown dark energy (or cosmologi cal constant ) whic h 



acts as a repu l sive fo rce (for reviews see iBahcall et al 
Frieman et al.ld2008l) ). 



1999), 



the first maximum of the power spectrum depends on the total 



Independently, the matter density parameter has been deter- 
mined from clustering of galaxies in the 2-degree Field Red- 



Supernova Cosmology Project 
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Figure 17: Combined constraints to cosmological densities Qa 
and Q.m, using supernovae, CMB and cluster abundance data. 
The flat Universe with Qa + Qm = 1 is shown with solid line 



(IKnopet al. . 2003) 



shift Survey and the Sloan Digital Sky Survey. The most ac- 
curate estimates of cosmological parameters are obtained us- 
ing a combined analysis of the 2dFGRS, SPSS and the WMAP 
data ( Spergel et al. , 2003 : Tegmarketal. , 2004 ; Sanchez et all 
2006). According to these studies the matter density parameter 
is Q m = 0.27 + 0.02, not far from th e value O m = 0.3, sug- 
gested by Ostriker & Steinhardtl ( 1995 ) as a concordant model. 
The combined method yields for the Hubble constant a value 
h = 0.71 + 0.02 independent of other direct methods. From 
the same dataset authors get for the density of baryonic matter, 
Q.b - 0.041 + 0.002. Comparing both density estimates we get 
for the dark matter density Qdm = Qm -&b = 0.23, and the dark 
energy density Qa = 0.73. These parameters imply that the age 
of the Universe is 13.7 + 0.2 Gigayears. 

6.2 The role of dark energy in the evolution of 
the Universe 

Studies of the Hubble flow in nearby space, using observations 
of type la supernovae with the Hubble Space Telescope (HST), 
were carried out by several groups. The major goal of the study 
was to determine the value of the Hubble constant. As a by- 
product also the smoothness of the Hubble flow was investi- 
gated. In this project supernovae were found up to the red shift 
(expa nsion speed) 20 000 km s _1 . This project ( Sandage et all 
2006) confirmed earlier results that the Hubble flow is very quiet 
over a range of scales from our Local Supercluster to the most 
distant objects observed. This smoothness in spite of the in- 



homogeneous local mass distribution requires a special agent. 
Vacuu m energy as the solutio n has been proposed by several au- 
thors dBarvshev et all d200ll) and others). Sandage emphasises 
that no viable alternative to vacuum energy is known at present, 
thus the quietness of the Hubble flow gives strong support for 
the existence of vacuum energy. 

7 Conclusions 

The discoveries of dark matter and t he cosmic web are two 
stage s of a typical scientific revolution dKuhnll 19701: Tremaine, 
1987). As often in a paradigm shift, there was no single discov- 
ery, new concepts were developed step-by-step. 

First of all, actually there are two dark matter problems - the 
local dark matter close to the plane of our Galaxy, and the global 
dark matter surrounding galaxies and clusters of galaxies. Dark 
matter in the Galactic disk is baryonic (faint stars or jupiters), 
since a collection of matter close to the galactic plane is possible, 
if it has formed by contraction of pre-stellar matter towards the 
plane and dissipation of the extra energy, that has conserved the 
flat shape of the population. The amount of local dark matter 
is low, it depends on the mass boundary between luminous stars 
and faint invisible stars. 

The global dark matter is the dominating mass component 
in the Universe; it is concentrated in galaxies, clusters and su- 
perclusters, and populates also cosmic voids. Global dark mat- 
ter must be non-baryonic, its density fluctuations start to grow 
much earlier than perturbations in the baryonic matter, and have 
at the recombination epoch the amplitude large enough to form 
all structures seen in the Universe. Initially neutrinos were sug- 
gested as particles of dark matter (hot dark matter), but presently 
some other weakly interacting massive particles, called cold 
dark matter, are preferred. 

Recently direct observational evidence was found for the 
presence of Dark (or vacuum) Energy. New data suggest that 
the total matter/energy density of the Universe is equal to the 
critical cosmological density, the density of baryonic matter is 
about 0.04 of the critical density, the density of dark matter is 
about 0.23 of the critical density, and the rest is dark energy. 

A number of current and future astronomical experiments 
have the aim to get additional data on the structure and evolu- 
tion of the Universe and the nature and properties of dark mat- 
ter and dark energy. Two astronomical space observatories are 
planned to be launched in 2008: the Planck CMB mission and 
the Herschel 3.5-m infrared telescope. The main goal of the 
Planck mission is to measure the CMB radiation with a precision 
and sensitivity about ten times higher than those of the WMAP 
satellite. This allows to estimate the values of the cosmological 
parameters with a very high accuracy. The Herschel telescope 
covers the spectral range from the far-infrared to sub-millimeter 
wavelengths and allows to study very distant redshifted objects, 
i.e young galaxies and clusters. 

Very distant galaxies are the target of the joint project 
GOODS - The Great Observatories Origins Deep Survey. Ob- 
servations are made at different wavelengths with various tele- 
scopes: the Hubble Space Telescope, the Chandra X-ray tele- 
scope, the Spitzer infrared space telescope, and by great ground- 
based telescopes (the 10-m Keck telescope in Hawaii, the 8.2- 
m ESO VLT-telescopes in Chile). Distant cluster survey is in 
progress in ESO. 

NASA-DOE have approved the mission DESTINY - the 



Dark Energy Space Telescope. Its goal is to detect and obtain 
precision photometry, light-curves and redshifts of over 2000 
type la supernovae over the redshift range 0.5 < z < 1.7 to con- 
strain the nature of dark energy. 

The largest so far planned space telescope is The James 
Webb Space Telescope (JWSP) - a 6.5-m infrared optimized 
telescope, scheduled for launch in 2013. The main goal is to 
observe first galaxies that formed in the early Universe. 

To investigate the detailed structure of our own Galaxy the 
space mission GAIA will be launched in 201 1. It will measure 
positions, proper motions, distances and photometric data for 1 
billion stars, repeatedly. Its main goal is to clarify the origin and 
evolution of our Galaxy and to probe the distribution of dark 
matter within the Galaxy. 

The story of the dark matter and dark energy is not over yet 
- we still do not know of what non-baryonic particles the dark 
matter is made of, and the nature of dark energy is also unknown. 
Both problems are a challenge for physics. So far the direct in- 
formation of both dark components of the Universe comes solely 
from astronomical observations. 
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